A comparison of sexual and asexual replication strategies in a simplified model based 

on the yeast life cycle 



Emmanuel TannenbaunJll 

Department of Chemistry, Ben-Gurion University of the Negev, Be'er-Sheva, Israel 

This paper develops simplified mathematical models describing the mutation-selection balance 
for the asexual and sexual replication pathways in Saccharomyces cerevisiae, or Baker's yeast. The 
simplified models are based on the single-fitness-peak approximation in quasispecies theory. We 
assume diploid genomes consisting of two chromosomes, and we assume that each chromosome 
is functional if and only if its base sequence is identical to some master sequence. The growth 
and replication of the yeast cells is modeled as a first-order process, with first-order growth rate 
constants that are determined by whether a given genome consists of zero, one, or two functional 
chromosomes. In the asexual pathway, we assume that a given diploid cell divides into two diploids. 
For the sake of generality, our model allows for recombination. In the sexual pathway, we assume 
that a given diploid cell divides into two diploids, each of which then divide into two haploids. The 
resulting four haploids enter a haploid pool, where they grow and replicate until they meet another 
haploid with which to fuse. In the sexual pathway, we consider two mating strategies: (1) A selective 
strategy, where only haploids with functional chromosomes can fuse with one another; (2) A random 
strategy, where haploids randomly fuse with one another. When the cost for sex is low, we find that 
the selective mating strategy leads to the highest mean fitness of the population, when compared 
to all of the other strategies. We also show that, at low to intermediate replication fidelities, sexual 
replication with random mating has a higher mean fitness than asexual replication, as long as the 
cost for sex is low. If the fitness penalty for having a defective chromosome is sufficiently high and 
the cost for sex sufficiently low, then at low to intermediate mutation rates the random mating 
strategy has a mean fitness that is a factor of \/2 larger than the asexual mean fitness. This is 
consistent with previous work suggesting that sexual replication is advantageous at high population 
densities, low replication rates, and intermediate replication fidelities. The results of this paper 
also suggest that S. cerevisiae switches from asexual to sexual replication when stressed, because 
stressful growth conditions provide an opportunity for the yeast to clear out deleterious mutations 
from their genomes. That being said, our model does not contradict theories for the evolution of 
sex that argue that sex evolved because it allows a population to more easily adapt to changing 
conditions. 
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I. INTRODUCTION 

The evolution and maintenance of sexual replication 
is one of the central problems in evolutionary biology 
[TJ m I3J m IS]. While sex is observed to be the pre- 
ferred, and often the only, replication strategy in the 
more complex, multicellular organisms, various forms of 
sexual replication are known to occur in single-celled or- 
ganisms as well. 

The ubiquity of sexual replication suggests that the 
replication strategy provides a selective advantage over 
asexual replication. The two basic theories for the se- 
lective advantage for sex are that sex provides a mech- 
anism for removing deleterious genes from a population 
[H m El m in] , and that sex allows for faster adaptation in 
dynamic environments IH ITJ [Sj • While these two ba- 
sic theories are often presented as mutually contradictory, 
this is not necessarily the case: That is, sex may provide 
a mechanism for removing deleterious genes from a popu- 
lation, and this same mechanism may also be responsible 
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for allowing faster adaptation in dynamic environments. 

The various theories for the existence of sex are in- 
complete, in that they provide an explanation for the 
selective advantage of the strategy, but do not explain 
why some organisms are essentially asexual, with some 
ability for recombination, while other organisms replicate 
exclusively sexually. 

In a recent set of papers [TTl [121 [131 [14] , Tannenbaum, 
Fontanari, and Lee argued that sex is generally favored in 
more slowly replicating organisms, since in this case the 
time (and energy) costs associated with finding a recom- 
bination partner are small compared to the characteristic 
growth time of the organism. For more quickly replicat- 
ing organisms, the time costs associated with sex become 
sufficiently large so that the asexual strategy becomes ad- 
vantageous. 

The results of the models presented by Tannenbaum, 
Fontanari, and Lee are broadly consistent with observa- 
tion. However, these models have the drawback that they 
assume highly simplified sexual replication pathways that 
do not exactly correspond to the sexual replication path- 
ways in actual organisms. These pathways were con- 
sidered because they were analytically solvable, and yet 
deemed a sufficiently good approximation of actual path- 
ways to yield biologically relevant results. Nevertheless, 
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a proper understanding of the evolutionary basis for sex 
will only be achieved when the asexual and sexual repli- 
cation pathways of actual organisms axe analyzed. 

In this paper, we seek to develop and analyze more re- 
alistic models comparing asexual and sexual replication 
strategies in unicellular organisms. In this vein, we are 
interested in constructing models based on the asexual 
and sexual rephcation pathways in Saccharomyces cere- 
visiae, or Baker's yeast. 

S. cerevisiae is a model organism that is used to investi- 
gate numerous fundamental problems in biology. Among 
these is the evolutionary basis for sexual replication. The 
reason for this is that S. cerevisiae, when stressed, en- 
gages in a form of sexual replication that is essentially 
the unicellular analogue of the gamete-based sexual repli- 
cation pathway in more complex organisms. Because S. 
cerevisiae is capable of both asexual and sexual repli- 
cation, understanding the reason for this organism to 
choose one replication strategy over another in a given 
situation will provide a clearer picture of the advantages 
and disadvantages of sexual replication over asexual repli- 
cation. 

This paper is organized as follows: In the following 
section (Section II) , we develop and analyze a model de- 
scribing the evolutionary dynamics of a population of 
S. cerevisiae replicating asexually. In Section III, we de- 
velop and analyze two models describing the evolutionary 
dynamics of a population of S. cerevisiae replicating sex- 
ually. The first model assumes a selective mating strat- 
egy, whereby only viable haploids may recombine with 
one another, while the second model assumes a random 
mating strategy. We find that when the cost for sex, 
as measured by the characteristic haploid fusion time, is 
negligible, then the selective mating strategy is advanta- 
geous over asexual replication for all replication fidelities, 
while the random mating strategy is only guaranteed to 
be advantageous over asexual replication at low to inter- 
mediate replication fidelities. In Section IV, we consider 
the case of self-fertilization, and show that, while such a 
strategy can be advantageous over pure asexual replica- 
tion, the sexual replication strategies of Section III are 
still advantageous at low to intermediate replication fi- 
delities, as long as the cost for sex is negligible. Finally, 
in Section V, we summarize the main results of this pa- 
per and describe plans for future research. In particular, 
we discuss how the results in this paper provide a plau- 
sible explanation for the benefits of sexual replication as 
a stress response in S. cerevisiae. 

It should be noted that the models we consider in this 
paper are highly simplified, in that we assume genomes 
consisting of only two chromosomes, and a fitness land- 
scape that is analogous to the single-fitness-peak land- 
scape for single-stranded genomes. Therefore, while 
these models are not intended to be quantitative at this 
stage, they are the two-chromosomed, single-fitness-peak 
analogues of the actual replication dynamics of S. cere- 
visiae, and arc therefore more realistic than previous 
models of unicellular sexual replication. 



II. ASEXUAL REPLICATION 
A. Definitions 

We begin our analysis by considering the evolutionary 
dynamics of a unicellular population replicating asex- 
ually. For simplicity, we assume that the genome of 
each cell consists of two chromosomes, and that, by 
analogy with the single-fitness-peak landscape in single- 
chromosomed models, a given chromosome is functional 
if and only if it is equal to some master sequence, denoted 
o-Q. It follows that a given genome consists of either zero, 
one, or two functional chromosomes. 

We assume that replication is characterized by first- 
order kinetics, and that the first-order growth rate con- 
stant depends on the number of functional chromosomes 
in the genome. We let n^y, k^u, and k„„ denote the first- 
order growth rate constants of cells with genomes consist- 
ing of two, one, and zero functional chromosomes, respec- 
tively. Here, "v" stands for "viable", and "u" stands for 
"unviable" . Naturally, we assume that i^yu ^ ^vu ^ 
Furthermore, for the remainder of this paper, we will as- 
sume that Kuu = 0, and we will also define a = Kvu/k.vv, 
so that a is the fitness penalty associated with having a 
defective chromosome. 

We also let n^y, riyu, and nyy denote the number of or- 
ganisms with two, one, and zero functional chromosomes 
respectively. We may then define the population frac- 
tions Xyy, Xyy, aud X yy aS X yy TiyyjTl, Xyy Tlyy^Tl, 

Xuu = nuu/n, where n = riyy + n„„ -|- n„„ is the total 
population. 

When a cell divides into two new diploids, the two 
parent chromosomes first replicate, and then distribute 
themselves evenly between the two cells. We assume that 
replication is not error-free, but that a given chromosome 
has a probability p of being replicated correctly. If we ne- 
glect backmutations, then a non-functional chromosome 
cannot produce a functional daughter. Therefore, a "v" 
chromosome produces a "v" daughter with probability 
p, a "u" daughter with probability 1 — p, while a "u" 
chromosome produces a "u" daughter with probability 
1. 

Finally, to allow for mitotic recombination within the 
limits of this model, we do not assume that the two 
chromosomes of a given parent-daughter pair segregate 
into distinct cells. Rather, if we tag one of the par- 
ent chromosomes with a "1", the other parent with a 
"2", the daughter of the "1" chromosome with a "3", 
and the daughter of the "2" chromosome with a "4", 
then the "1" chromosome has a certain probability of co- 
segregating with any of the other chromosomes. Thus, we 
let ri2 denote the probability that the two parent chro- 
mosomes co-segregate into one of the daughter cells, so 
that the two daughter chromosomes co-segregate into the 
other daughter cell; we let ri3 denote the probability that 
parent chromosome "1" co-segregates with its daughter, 
so that parent chromosome "2" co-segregates with its 
daughter; and we let ri4 denote the probability that par- 
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FIG. 1: (Color online) Illustration of the various replication 
pathways and their associated probabilities, as a function of 
the recombination probabilities ri2, ria, and ru. 



ent chromosome "1" co-segregates with the daughter of 
parent chromosome "2". Note that ri2 + ri3 + ri4 = 1. 



B. Population genetics equations 

With the recombination probabihties ri2, ri3, ru and 
repHcation fidehty p in hand, we may compute the prob- 
abihties of the various rephcation pathways, as shown in 
Figure 1. This leads to the following system of differen- 
tial equations governing the time evolution of the various 
population fractions: 



= [Kyy{2p + ri2(l - p)2 - 1) - R{t)]3 



dt 

[k„„p(1 - 2ri3) - K(t = oo)]xy 



+2KyvXvv{i -p)[p+ (1 - ri2){l -p)] 

+ KvvXyvri2{l - pY 



(1) 



We have introduced an additional parameter, R{t), 
which is the mean fitness of the population and is de- 
fined as K,{t) = {\/n){dn/dt) = Ky^Xw + K„„a;„„. The 
mean fitness is the per capita growth rate of the popula- 
tion, and therefore measures the first-order growth rate 
constant of the population as a whole. 

In the group selection approach that will be adopted 
in this paper, the central object of study is the steady- 
state mean fitness of a population. The reason for this 
is that if two or more populations are growing exponen- 
tially, the population will the largest mean fitness will 
drive the others to extinction. Therefore, for a given set 
of parameters, the replication strategy with the largest 
steady-state mean fitness is the advantageous one. 



C. Mutation-selection balance and the steady-state 
mean fitness 

It is a standard result from population genetics that 
the above system of differential equations will converge 
to a steady-state, corresponding to what is known as a 
mutation- selection balance. Therefore, to determine the 
mutation-selection balance, we set the left-hand-side of 
the above system of equations to 0, giving that, 

= [nvvi^p + ri2(l - p)^ - 1) - K.{t = co)]xyy + Kyurispxvu 
= [uvuPi^ - 2ri3) - K{t = oo)]xvu 
+2KvyX.uy{l -p)[p+ (1 - ri2)(l ~p)] 

= -R{t = CO)XuU + KyyXyui^ ^ ^ ^"13)^] 

{l-pf (2) 
Solving the first equation for Xw in terms of Xyu gives, 

" R{t^oo)~Kyy{2p + r,2il-pr-l) ^ ' 
Substituting into the second equation gives, 

= Xyu[KyuP{l - 2ri3) - R{t ^ oo) 
+2(l-p)b+(l-ri2)(l-p)] X 

R{t ^ oo) ~ Kvv{2p + ri2(l — p)'^ — 1) 

and so, we either have that Xyu = Xyy = ^ Xuu = 1, 
so that R{t = oo) = 0, or R{t — oo) is the solution to the 
quadratic, 

^ K(t = oo)2 - [K„„(2p + ri2(l - 1) 
+Kyup{l - 2ri3)]R{t = oo) 

+ KyyKyyP[2{l " ri3)p + TjaCl - ^ l] (S) 

Dividing both sides by Kyy, we obtain that the nor- 
malized steady-state mean fitness, (f) = R(t = oo)/Kyy, 
satisfies the quadratic, 

= [2p + ri2(l-p)2-l + ap(l-2ri3)]0 

+ap[2(l-ri3)p-Hri2(l-p)2-l] (6) 

as long 3jS X ufj and Xyu are not both 0, which is equivalent 
to the statement that > 0. 

The quadratic in Eq. (6) has two solutions. To deter- 
mine which of these solutions corresponds to the mean 
fitness of the population, we examine the two solutions 
for both p = and p ~ I. When p — 0, Eq. (6) becomes. 







'+{l-ri2)] 



(7) 



which admits the solutions = and (p = ri2 — 1 < 0. 
The physical solution is (j) — 0, since a negative mean 
fitness is impossible for this model (this model does not 
include death). 

When p = I, Eq. (6) becomes. 



0= [0-a(l-2ri3)][0-l] 



(8) 



which admits the solutions (j) = a(l — 2ri3) < 1, and 
(j) = 1. The physical solution is (/> = 1, since a population 
should consist entirely of the wild-type when replication 
is error-free. 

Note that for both p = and p = I, the value of (p is 
the larger of the two roots of the quadratic polynomial 
in Eq. (6). By continuity, we expect that the value of (j) 
for all p e [0, 1] is given by the larger of the two roots of 
the quadratic polynomial in Eq. (6). 

We may also evaluate (j> for a = and a = 1. When 
q; = 0, we obtain that (p — max{2p -I- 7-12(1 — p)^ — 1, 0}, 
while when a = 1 we have that (j) is the solution to, 

0= [0-p][0-(2(l-ri3)p + ri2(l-p)2-l)] (9) 

which admits the solutions (j) ~ p and (f) — 2(1 — ri3)p + 
^■12(1 - P)^ - 1- Because 2(1 - rig)^ -I- ri2(l - p)^ - I < 
1 = < p, we have that (f> — p for a = 1 . 



2p+(l-p)2- 
When ri3 = 



0, we obtain that 6 is the solution to. 



0= [0-ap][0-(2p+ri2(l-p)' 



(10) 



so that (f> — max{2p -I- 7-12(1 — p)^ — 1, ap}. 

Finally, when a,p, ri3 € (0,1), we can show that (j) G 
{ap,p). To do this, note first that the polynomial of 
Eq. (6) goes to 00 as (/i —00. When <j) = ap, the 
polynomial evaluates to — 2ri3a(l — a)p^ < 0, and so, by 
the Intermediate Value Theorem, the polynomial of Eq. 
(6) has a root in (—00, ap). 

When 4> — P; the polynomial evaluates to p(l — p)[l -I- 
a — (1 — Q!)ri2(l— p)] > 0, and so, again by the Intermedi- 
ate Value Theorem, the polynomial has a root in {ap,p). 
Because the polynomial of Eq. (6) has only two roots, 
and because the larger of the two roots is the value of (j), 
we have that 4> G i^ip^p) whenever a,p, ri3 G (0, 1). 



D. Limiting behaviors of 

We may determine (d0/rfp)p=o,i by differentiating 
both sides of Eq. (6) and substituting p = and p = I, 
respectively. For p = 0, we obtain, after some manipula- 
tion. 



d(j) 
^dp>^=' 



(11) 



while when p = 1 we obtain. 



' dp 



)p- 



^ 1 - q;(1 - ri3) 
'l-a(l-2ri3) 



(12) 



Note then that when p is close to 0, we have that 
ap, independently of the values of ri2, ri3, and ri4. 



III. SEXUAL REPLICATION 
A. Definitions 

The sexual pathway in yeast involves the division of 
a diploid into two daughter diploid cells, each of which 
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FIG. 2: Illustration of the sexual replication pathway in S. 



cerevisiae. 



divides again to produce a total of four haploids. The 
four haploids then enter a haploid pool, where they con- 
tinue replicating until they encounter another haploid 
with which they fuse to form a diploid that may grow 
and then repeat the cycle. This process is illustrated in 
Figure 2. 

The haploid cells of S. cerevisiae belong to one of two 
distinct mating types, termed a and a. A given haploid 
can only fuse with another haploid of a different mating 
type, so that only a — a matings are possible. However, 
when a haploid of a given mating type divides, the daugh- 
ter can sometimes spontaneously switch mating types. 

To develop a model for sexual replication in yeast, 
we need to define a few additional quantities. To this 
end, we let k„ and k„ denote the first-order growth rate 
constants of viable and unviable haploid yeast, respec- 
tively. As with diploid yeast, we assume that k„ = 0. 
We let n.u^a, n^,a, 7T.«,a, n^^a denote the number of vi- 
able a, viable a, unviable a, and unviable a haploids 
in the population. We also define n = n^y + Uyu + 
riuu + {l/2)[ny^a + riy^a + "«,a + n«,a]- Finally, we de- 
fine the population fractions Xyy = nyy/n, Xyu = Uyy/n, 

•^uu ^uul^i -^v^a ( l/2)7T'i;^(j/71, XyQ, {\ j^^Tiy / Tl, 
Xu,a = {l/2)ny^a/n, Xu,a = {l/2)nu,a/n. 

We assume that haploid fusion may be modeled as a 
second-order rate process with rate constant 7. We also 
assume that, as the population grows, the system volume 
expands to maintain a constant population density p, de- 
fined as p = n/V. Finally, we let Pm denote the prob- 
ability that when a haploid cell replicates, the daughter 
will be of the opposite mating type from the parent (be- 
cause S. cerevisiae replicates by budding, it is possible to 
canonically define a single daughter cell). 

We consider two distinct mating strategies: (1) Selec- 
tive mating, where we assume that only viable haploids 
may fuse with one another. (2) Random mating, where 
all haploids may participate in the mating process, and 
where haploids do not exhibit any kind of preference for 
a certain genome type. 

The reasoning behind considering a selective mating 
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strategy is two-fold: On the one hand, it may be argued 
that haploids with unviable genomes are simply physi- 
cally unable to participate in the mating process. On the 
other hand, it may be that the viable haploids have some 
mechanism for discerning the genome of another haploid 
(say by reading what are known as indicator traits) |15j . 
and essentially choose to only fuse with another viable 
haploid. 

The reasoning behind considering a random mating 
strategy is that in practice, it is very difficult to deter- 
mine another organism's genome. Indicator traits (i.e. 
the phenotype) are not in one-to-one correspondence 
with the genome (i.e. the genotype). Furthermore, the 
actual process of determining an organism's fitness takes 
time and consumes energy, all of which incur additional 
fitness penalties on the sexual replication strategy. These 
two factors can make the selective mating strategy highly 
impractical if not impossible. The random mating strat- 
egy does not suffer from the difficulties of the selective 
mating strategy, and so it is important to consider its 
selective advantage, if any, over asexual replication, if 
we are to draw biologically relevant conclusions from our 
two-chromosome model. 

The difficulty and the fitness penalty associated with 
reading an organismal genome exactly on the one hand, 
and the benefit of fusing with a fit haploid, likely means 
that actual mating strategies in real organisms lie some- 
where between the two extremes being considered here. 
Nevertheless, with the highly simplified models we are 
considering, the selective and the random mating strate- 
gies are the most natural ones to initially study. 



B. Selective mating 

The population genetics equations for the selective 
mating strategy may be readily shown to be, 



state equations, 

= -(1 + cb)x,, + 

= —(j)Xv + Xvvi'^ +p) H —pxv — -^xl 

= ~(t)Xu + Xvy{l - p) + —{1 - p)Xv (14) 

We will now consider the solution to these equations 
in the limit where ^p/nyv oo. In this regime, the 
contact rate between the haploids is so large that Xy ^ 
at steady-state. Adding the first two of the previous set 
of equations, and setting Xy — 0, gives, 

= Xyy{p-(f)) (15) 

so that (f> — p a.t steady-state. 

Clearly, then, the selective mating strategy has a 
higher mean fitness than the asexual replication strat- 
egy, assuming that ^p/nyv oo, or equivalently, that 

Kvv/ilP) 0. 

Note that Kyy/^jp) measures the ratio of the charac- 
teristic fusion time, given by l/(7p), to the characteristic 
growth time of the cell, given by 1 When this ratio is 
small, that is, when the fraction of time a haploid spends 
looking for a mate with which to rccombine is small com- 
pared to the length of the cell cycle, then the steady-state 
mean fitness of the population is larger than that of the 
corresponding asexually replicating population. 

The quantity Kyy/^jp) is referred to as the cost for 
sex in this model, for when it is small then there is little 
fitness penalty associated with sexual replication, while 
when Kyy/ (7p) is large, then jp/ Kyy is small, and so hap- 
loid fusion is slow, leading to a possibly significant fitness 
penalty associated with the sexual replication strategy. 

What we have established is that, when the cost for 
sex is negligible, then a selective mating strategy leads 
to a higher steady-state mean fitness than the asexual 
strategy. 



dt 

dXyi 

~dt 
~dt 



K{t)]Xyy + ^PXI 
- K{t)]Xyu 



/s^ Xii 



1 



'K(^t)Xy -\- KyyXyy{l -\- p) -\- ^KvuXyu{^ ~^ P) 



+ KyXyP - JPX^ 

dx„, 1 

-Hl{t)Xn -|- KyyTLyyiX p) + „ Ti^, ^ ( 3 p) 



dt 

+ KyXy{l -P) 



(13) 



Note that since there is no production of vu or uu 
genomes, we may assume that = x^^ = at steady- 
state. This gives, after normalizing by Kyy, the steady- 



C. Random mating 

For random mating, the population genetics equations 
may be readily shown to be, 



dt 

d/X 

~dt 

dxi,i 



-{Kyy + K(t))Xyy + J pxl 
~{l^VU + ^{t))Xyii -\- 1^ pXyXy 

-K{t)Xuu + IPxl 



dt 

dXy 1 

— ^ = K{t)Xy + (1 + P^KyyXyy -j" -^KyuXyuiX "I" P) 
-\-KyXyP 'ypXy^Xy -\- X y) 

dXy, 1 

^(^)^U (-^ P)^VvXyy -\- -zKyiiXyy^ij^ p) 



dt ^ ' ' — - 2 

+ KyXy{l -P) - JpXu{Xy + X y) 



(16) 
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If wc set the left-hand-sides of the above equations to 
and normalize by Kyy, we obtain that the steady-state 
equations are, 

= -(l + ,^).T„, + (^)x2 

= -{a + 4>)xyu + 2{—)x^Xu 

Q = -(pXuu + { — )xl 

1 n 

Z f^vv 

1 K 

= -^Xu + (1 -p)a;^^ + ^o;(3 H -p) 

( ^Xy,(^Xy -\- Xiij 

Kyy 



Now, as ^pji^yy oo, we obtain that Xy^x^ 
so, adding the last two equations, we obtain, 

{:i^){x, + xuf = 2cp 
Defining Xy = Xv/{xy + Xu), we then have, 



(1 + P){Xyy + -aXyu) = 2(/)i„ 



(17) 
0, and 

(18) 
(19) 



Plugging this expression into the first two equations, 
we obtain, 

ct>il + 4>)xyy = {l+pf^^±^ 

(j){a + (t)){(f)- Xyy) = + P)4>{4> + Xyy) 



(20) 



The first equation may be used to replace (</> + Xyy)"^ 
in the second equation by a term that is proportional to 
Xyy. The second equation may then be used for solve for 
Xyy. The result is. 



Xyy 



(f> — ap 



1 - 2a) 



ap 



(21) 



This expression may then be plugged back in to the first 
equation, which gives, after some manipulation, that, 

+ 1] [(/. - ap] [(/>(! - 2a) + ap] - ^ (1 + p)2(l - a)2,/,2 ^ 

(22) 

This expression defines the steady-state mean fitness of 
the sexual population with random mating, when the cost 
for sex is negligible. 

When (j> = ap, the left-hand-side evaluates to 
-(1/2)(1 -|-p)2(l - a)2aV < 0, while when (/) = p, the 
left-hand-side evaluates to (l/2)(l-p^)p^(l-a)^ > 0, so 
that (j> has a solution in [ap,p], which is the normalized 
steady-state mean fitness of the population. 



When p = we obtain that (j) = 0, while when p = 1 
we obtain that (p = 1. When a = we obtain that 
(f) — max{(l/2)(l + p)2 — 1,0}, while when a = 1 we 
obtain that (j) = p. 

Differentiating both sides of Eq. (22) twice, and sub- 
stituting p = 0, gives, after some manipulation. 







2a^ (23) 



and so. 



,d(t) 



a+lTV2 



(g)..o = a(2±V2)-^,^^^_^ 



(24) 



Now, 0^ + 20.-1 = {a + l + s/2){a + l-s/2), so that. 



# /- V2±l 

[-r)p=o = v2a = 

dp'^ a+l±V2 



(25) 



Since we expect (() to be positive for p > 0, we take the 
+ solution, giving. 



2 + V2 
+ 1 + V2 



(26) 



Note that (2 + \/2)/(a+l + \/2) decreases from ^2 > 1 
to 1 as a increases from to 1. Therefore, when p is close 
to 0, (6 « [(2 + V2)/{a + 1 + V2)]a:p > ap for a < 1, so 
that the mean fitness of the sexual population is greater 
than that of the asexual population when the replication 
fidelity is at a low to intermediate value. 

Because ([(2 + y^) /{a+l + V2)] — l)ap increases with 
p, the fitness difference between sexual and asexual repli- 
cation initially increases as p increases away from 0. If 
the cost for sex is non-negligible, sexual replication with 
random mating will only have a selective advantage over 
asexual replication if this initial fitness difference is suf- 
ficiently large. This implies that sexual replication with 
random mating will only have a higher mean fitness than 
asexual replication at intermediate replication fidelities. 

Differentiating both sides of Eq. (22) and substituting 
p = 1 gives, after some manipulation. 



(27) 



As a final note for this subsection, we shall prove that 
the sexual mean fitness with random mating is greater 
than the asexual mean fitness when ri2 = t"i4 = 1/2, 
ri3 = 0, a,p G (0, 1), and there is no cost for sex. 

We begin by defining (f)a to be the normalized mean 
fitness for the asexual population, and 0,.,, to be the nor- 
malized mean fitness for the sexual population with a 
random mating strategy. 

First of all, when p = 0, then (j)rs = (pa = 0, and 
when p = 1, then cprs = = 1. When a = 0, then 
<Aa = (prs = niax{2p -|- ri2(l — p)^ — 1,0}, while when 
a = 1, then = = p. So, in what follows, we will 
assume that a,p € (0, 1). 
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We have that {d<j)a/dp)p^i = 2[1 - a{l - ri3)]/[l - 
a(l — 2ri3)], so that, when ri3 > 0, {dcpa/ dp)p=i <2 — 
{d4)rs/dp)p=i. This impHcs that (pa > (f'rs for p sufh- 
ciently close, but not equal, to 1. Since we also have that 
(j}rs > 4'a for p sufficiently close, but not equal, to 0, this 
implies that (pa is initially larger than (prs for p close to f , 
but then must equal at some p e (0, f ) and eventually 
remains below (p^s for p sufficiently close to 0. 

Now, when ri3 — 0, we fcnow that (pa = max{ap, 2p + 
''12(1 — ~ !}• I3ecause 4>rs > cup for a,p e (0,f), 
we need only show that 4>rs > 2p + ri2(f — p)^ — 1 for 
ri2 = 1/2. 

Suppose that 0,.^ = 2p+ (f/2)(f -p)^ - f = (f/2)(f + 
p)^ — 1 for some p e (0, f ). Substituting into Eq. (22) we 
obtain, after some manipulation, that ^j,^ = p. But then 
we have that (1/2) (f + p)^ — 1— p^p=l, which is a 
contradiction, since p G (0, 1) by assumption. 

Therefore, (prs + 2p+ (f /2)(f -p)^ - f for a,p e (0,1), 
and so <prs > 4>a whenever a,p £ (0,1), ^13 — 0, and 
ri2 = ri4 = 1/2, as we wished to show. 



IV. SELF-FERTILIZATION 

A final replication strategy that we will consider is a 
form of self-fertilization. Here, the four haploid cells pro- 
duced from a given diploid simply fuse with one another. 
We consider this replication mechanism independently of 
whether or not it occurs in S. cerevisiae. The reason for 
this is that we want to examine whether the selective ad- 
vantage obtained in our sexual replication models could 
have been obtained without needing to fuse with hap- 
loids from a distinct diploid cell. If it turns out that 
self-fertilization yields a mean fitness at least as high as 
the one obtained from our sexual models, then the mod- 
els we have considered thus far do not provide a clear 
picture of the selective advantage of sexual replication 
over asexual replication. 

The random mating strategy is equivalent to the asex- 
ual replication strategy where ri2 — ri3 — rn = 1/3. 
Since we already know that the sexual strategy with ran- 
dom mating will outcompete any asexual strategy at suf- 
ficiently low replication fidelities, it follows that the sex- 
ual strategy with random mating will outcompete self- 
fertilization with random mating at sufficiently low repli- 
cation fidelities. 

We now turn to the case of selective mating. Here, 
we assume that only two "v" haploids can fuse with one 
another. However, if out of a group of four haploids pro- 
duced from a given diploid, there is only one "v" haploid 
present, we assume that that haploid is lost, and does 
not participate further in the replication process. 

The differential equations governing the time evolution 



of the population numbers n^v, riyu, and n^u is then, 

dflyv 2 , 

r; — KyyTlyyP -f- ^yu^vuP 

dt 

dlT'yu 

dt 

^ = (28) 

Note then that when the mutation-selection balance is 
reached, we have Uyu — n^u — 0, so that, in the long- 
term, we have n — Uyy, and dn/dt — Kyyp^n, which 
implies that (p — p^ < p ioi p G (0, 1). Therefore, self- 
fertilization with selective mating has a lower mean fit- 
ness than sexual replication with selective mating. 



V. DISCUSSION AND CONCLUSIONS 

The two key results of this paper are that, when the 
cost for sex is low, the sexual replication model with selec- 
tive mating has a higher steady-state mean fitness than 
the asexual replication model, and the sexual replica- 
tion model with random mating has a higher steady-state 
mean fitness than the asexual replication model at low 
to intermediate replication fidelities. Interestingly, for 
random mating, the steady-state mean fitness at low to 
intermediate replication fidelities goes from being equal 
to the asexual mean fitness when a = 1, to being a factor 
of \/2 larger than the asexual mean fitness as a decreases 
to 0. 

The conclusions of this paper provide a possible expla- 
nation for why S. cerevisiae engages in sexual replication 
as a stress response: When conditions are favorable, the 
replication rate, as measured by Hyy, is fairly large, and 
so the cost for sex, as measured by Kw/iip), is large as 
well, so that sex is the disadvantageous strategy. How- 
ever, when the yeast population is under stress, Kyy be- 
comes small, and so the cost for sex becomes small as 
well. If the mutation rate in yeast is not too low, then sex 
can become the advantageous strategy. This explanation 
is consistent with experimental work on yeast suggesting 
that sexual replication in yeast provides a mechanism for 
removing deleterious mutations from the yeast genome 

In a previous series of papers |lll I12j . Tannenbaum 
and Fontanari considered a different sexual replication 
strategy than the one being considered here. Tannen- 
baum and Fontanari assumed that a given diploid splits 
directly into two haploids, which enter a haploid pool, 
fuse, and then the resulting diploid divides asexually. For 
both the selective and random mating strategies, it was 
shown that, when the cost for sex is low and if ri3 — 0, 
sexual replication leads to a higher mean fitness than 
asexual replication. 

Although the replication models considered in [TTl [12] 
led to interesting results, they had two main problems: 
First of all, for the random mating strategy, the sexual 
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replication process becomes identical to asexual replica- 
tion when ri3 = 1. This limited the generality of the 
conclusions drawn regarding the selective advantage of 
sexual versus asexual replication in various parameter 
regimes. 

Second, although the models being considered in this 
paper are highly simplified, they are nevertheless based 
on a sexual replication pathway that is much closer to the 
actual sexual replication pathway in S. cerevisiae (within 
the constraints of the model). The replication models 
considered in [TTJ [T^], by contrast, were only loosely 
based on the sexual replication pathway in S. cerevisiae. 
As a result, we argue that the results in [TTJ [T^] could 
only be used to draw general conclusions about the selec- 
tive advantage for sexual replication. On the other hand, 
we argue that the results in this paper, although based 
on a highly simplified model, provide a starting point for 
understanding the selective advantage for sexual versus 
asexual replication in actual unicellular organisms, e.g. 
5*. cerevisiae. 

For future research, we would like to develop more re- 
alistic models involving multi-gened, multi-chromosomed 
genomes. Furthermore, this paper explored the selective 
advantage for sex in static environments. For future re- 
search, we would like to develop models describing com- 



petition between sexual and asexual replication in dy- 
namic environments, to test the theories for sex that are 
based on the idea that sex increases adaptability. 

In this vein, it should be noted that there has been 
experimental work on yeast suggesting that sex is advan- 
tageous because it allows yeast to adapt more quickly to 
changing environments [17!. As we have mentioned pre- 
viously, this view is not necessarily contradictory to the 
results of |16j , which argue that sex is a way of removing 
deleterious mutations. In our opinion, at intermediate 
mutation rates, and when the cost for sex is sufficiently 
low, there is an advantage for sex that holds in static 
environments and in an infinite population. This advan- 
tage, however, can be significantly enhanced in dynamic 
environments, and for intermediate population sizes 
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